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On the Continuum Limit of Topological Charge Density Distribution 
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The bulk distribution of the topological charge density, constructed via HP 1 cr-model embedding 
method, is investigated. We argue that the specific pattern of leading power corrections to gluon 
condensate hints on a particular UV divergent structure of HP 1 <r-model fields, which in turn im- 
plies the linear divergence of the corresponding topological density in the continuum limit. We show 
that under testable assumptions the topological charge is to be distributed within three-dimensional 
sign-coherent domains and conversely, the dimensionality of sign-coherent regions dictates the lead- 
ing divergence of the topological density. Confronting the proposed scenario with lattice data we 
present evidence for indeed peculiar divergence of the embedded fields. Then the UV behavior of 
the topological density is studied directly and is found to agree with our proposition. Finally, we in- 
troduce parameter-free method to investigate the dimensionality of relevant topological fluctuations 
and show that indeed topological charge sign-coherent regions are likely to be three-dimensional. 

PACS numbers: 11.15.-q, 11. 15. Ha, 12.38.Aw, 12.38.Gc 



I. INTRODUCTION 

Topology investigations had always been conspicuous 
topic for the lattice community and the recent advances 
indeed put it on the solid grounds. Both the topologi- 
cal charge and its density are now computable on ther- 
malized vacuum configurations and the results already 
obtained in pure Yang-Mills theories indicate that the 
conventional instanton based models are to be strongly 
modified. Here we basically mean the discovery of global 
topological charge sign-coherent regions and, what is 
even more important for us, the lower dimensionality of 
these regions 0, 0, 0, El IS E| > which are now believed 
to be three-dimensional. Note that the lower dimension- 
ality of physically relevant vacuum fluctuations should 
not come completely unexpected, it had been repeatedly 
discussed in the recent past (see, e.g., Refs. 0,111 El)- 

It is important that qualitatively the same picture 
of relevant topological fluctuations appeared recently 
within radicallydifferent approach introduced and devel- 
oped in Refs. (lCllllj . Without mentioning all the details 
and technicalities involved we note only that the HP 1 em- 
bedding method is the nearest to the classical ADHM in- 
vestigation of SU (2) gauge fields topology and essentially 
reconstructs the topology defining map S 4 — > HP 1 = S 4 , 
in terms of which both the topological charge and its den- 
sity obtain unambiguous and well defined meaning. We 
are not in the position to review all the results obtained 
in |ToL fTTj . however, it is important that HP 1 embedding 
allows to get rid of leading perturbative divergences in 
various observables. Moreover, it reproduces the topo- 
logical aspects of SU(2) Yang-Mills theory, is fairly com- 
patible with other topology investigation approaches and 
allows to calculate with amazing accuracy the gluon con- 
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densate, demonstrating firmly that its quadratic power 
correction does not vanish. It is crucial that the only 
non-trivial power dependence of HP 1 projected curva- 
ture upon UV cutoff is contained in the quadratic term, 
which must be encoded in the local structure of the topol- 
ogy defining map S 4 — > HP 1 . In section Hi Al we argue 
that the only possibility to reproduce the observed pat- 
tern of power corrections is to assume that the mapping 
S 4 — > HP 1 is highly asymmetric so that the correspond- 
ing Jacobian (which is essentially equivalent to the topo- 
logical density) diverges linearly with diminishing lattice 
spacing. In turn this divergence implies rather peculiar 
geometry of the relevant topological excitations. Namely, 
we show in section Hi Bl that the topological charge sign- 
coherent regions are to be three-dimensional domains em- 
bedded into original four-dimensional space. Here the 
consistency check is provided by the fluctuations of topo- 
logical charges associated with sign-coherent regions. 

Attempt to confront the above scenario with the lattice 
data brings out a wealth of both technical and theoret- 
ical issues, which are addressed in section IIIII In sec- 
tion llll Al we present numerical measurements which jus- 
tify our conclusion about the local structure of the map 
S 4 — > HP 1 . The direct measurement of the UV behavior 
of the topological density (sections IIII Bl IIII C|) necessi- 
tates the invention of new calculation algorithm, which 
turns out to be fast and rather accurate. It allows us 
to make statistically significant comparison of the theo- 
retical considerations with lattice data and confirms that 
the characteristic topological density is indeed divergent, 
but at most linearly, in the continuum limit. As far as 
the dimensionality of relevant fluctuations is concerned, 
we investigate it in section Till Dl and develop unambigu- 
ous, in fact, method of its determination. Essentially it is 
the specially crafted biased random walk model embed- 
ded into the ambient topological density environment [Tsj . 
We argue that the appearance of critical-like behavior of 
the model signifies the lower dimensional long-range or- 
der present in the topological density background. We 
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show that critical regime does occur and indicates that 
indeed the density is concentrated in highly extended 
submanifolds, the dimensionality of which is compati- 
ble with three. Finally we investigate the fluctuations 
of the topological charges associated with sign-coherent 
regions and argue that only three-dimensional domains 
with linearly divergent topological density are consistent 
with theoretical and experimental restrictions imposed 
on the structure of topological fluctuations. 

II. SCALING OF TOPOLOGICAL DENSITY 

A. Gluon Condensate, Leading Power Correction 
and Divergence of Topological Density 

The essence of HP 1 a- model embedding approach is 
the assignment of unique configuration of HP 1 cr-model 
fields \q x ) to every given SU(2) gauge background 
(until section IlIII we use continuum notations), where 
| q x ) is two-component, normalized (q x \q x ) = 1, quater- 
nionic vector (see Refs. [ldLLl| for details). The relevant 
configuration \q x ) is defined by the requirement that it 
provides the absolute minimum to the functional 

F(A,q) = J ATrfA, + (q\d^\q)) 2 (1) 

for given (fixed) gauge potentials A^. Note that gauge 
covariance is maintained exactly since the cr-model tar- 
get space HP 1 is the set of equivalence classes with re- 
spect to \q x ) ~ \q x )v x , v x € SU(2). The uniqueness 
of the minimum of and the factual absence of Gribov 
copies problem was discussed in length in Refs. 0, [n| 
and here we take it for granted. Therefore the embedded 
HP 1 cr-model fields I q x ) arc unique (although non-local) 
functions of the original potentials. The advantage of the 
construction is that the gauge fields topology becomes ex- 
plicit in terms of | q x ) variables. Indeed, the gauge invari- 
ant projectors | q x ) ( q x | provide the map of compactified 
physical space S 4 into the target space HP 1 , the degree 
of which is equal to the topological charge of the original 
gauge background. Furthermore, the local distinction of 
this map from trivial one is the uniquely defined measure 
of the topological charge density. One could say that the 
sole purpose of HP 1 embedding method is to locally re- 
construct the topology defining map S 4 — > HP 1 , which 
allows to get essentially all the topological aspects of the 
original background. 

Then it is natural to introduce the HP 1 projection 

A, - Af = -{q\d,\q), (2) 

which replaces A^ with its best possible approximation 
by cr-model induced potentials. The striking properties 
of the projected fields A^ p were investigated in details 
in Ref. ^lj. In particular, it was shown that A^ p ex- 
actly reproduce the most of non-perturbative aspects of 
the original SU(2) configurations, while containing no 



sign whatsoever of usual perturbative divergences. To 
the contrary, the kernel of the map J3J) was shown to 
correspond to pure perturbation theory with identically 
trivial topology and vanishing string tension. Without 
mentioning all the aspects and properties of HP 1 pro- 
jection @ let us note that it allows to calculate with 
inaccessible so far accuracy the gluon condensate and 
its leading power correction. Indeed, the spacing depen- 
dence of HP 1 projected gauge curvature was found 
to be astonishingly well described by 

<^(^) a > = ^ + ^<-G 2 >, (3) 

ol 2 = [61(3) MeV] 2 , (4) 
(^G 2 ) = 0.0271(10) GeV 4 , 

TT 

where a is the lattice spacing (1/a servers as the UV 
cutoff). It is crucial that Eq. © does not contain any 
sign of usual perturbative contribution of order 0(a~ 4 ). 
Note that the value of a>2 coefficient turns out to be 
unexpectedly small, nevertheless it fits nicely into the 
known bounds on the magnitude of the quadratic cor- 
rection term (see Ref. for recent review). While the 
actual numbers quoted in (@J) are not important for the 
present discussion, it is crucial that they both are defi- 
nitely non-zero and therefore the spacing dependence of 
(Tr(F^f) 2 ) includes only 0(a" 2 ) and O(a ) terms. 

To analyze the consequences of Eq. © let us consider 
point n £ 5 4 , in the neighborhood of which the map 
S 4 — > HP 1 is non-degenerate. In the vicinity of n and its 
image m £ HP 1 we introduce local coordinates x^ and 
// = 0, ...,3 (stereographic projection from points n 
and m correspondingly) such that y Pi (x = 0) = 0. Non- 
degeneracy means that det[dy/dx] ^ and hence the 
function y(x) is invertible. From the specific explicit form 
of HP 1 projected fields @ we conclude that 

Af{x) = -{q x i^i qx ) = j;wAr f (!/W), (5) 

where Jjf = dy^/dx" is the Jacobian matrix, A™ st {y) is 
the potential of classical BPST instanton solution with 
unit radius 

\ (y) = -(q\d»\q) = 2 (l + y 2 ) () 

and denotes quaternionic basis. Note that Eq. JSJ) is 
local and crucially depends upon the space-time varying 
matrix J^. Furthermore, Eq. (JSJ) relies heavily on the 
HP 1 projection and would not be valid for generic 
gauge potentials. The corresponding HP 1 projected cur- 
vature is similarly expressible in terms of instantonic 
field-strength 

pHP . jpi jp 2 rpinst (n\ 
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and therefore at point n € S we have 

\^{F^f cx (Tr 5 ) 2 - Tr 5 2 , (8) 

where we have introduced the metric g^ v — J^J\ and 
skipped inessential numerical factor. Therefore the study 
of the projected curvature ( \ Tr(F^ p ) 2 ) is equivalent to 
the local investigation of the metric g associated with 
the topology defining map «S 4 — ► HP 1 . In term of strictly 
positive eigenvalues A M of g, Eqs. 10 and (JSJ imply 

£ (^u) = A + B-a- 2 , (9) 

A ~ A QCD , i? ~ ^QCD ' 

where we have generically indicated the IR physical scale 
involved in Eq. J3J) by Aqcd and explicitly kept all pow- 
ers of lattice spacing. Evidently, Eq. imposes rather 
stringent restrictions on the distribution of eigenvalues 
A M and requires them to depend highly non-trivially upon 
the UV cutoff. However, Eq. Q is not sufficient to an- 
alyze this dependence in details. The relation, which 
provides an additional input and which is verified nu- 
merically in section Ull Al reads 

£(A M ) = a + f3-a- 2 , (10) 

a ~ A QCD > P ~ 1 • 

Taken at face value it indicates that generically all the 
eigenvalues are quadratically divergent in the limit a — > 
0. However, it turns out that the simplest ansatz 

(A )=/3-a- 2 , (Xi)=ai, i=l,2,3, (11) 

a i ~ A QCD > Hi a i = a > 

is capable not only to reproduce the observed pattern 
of power corrections |[?JJ, (|l(Jf> . but also passes stringent 
consistency check, which we describe next. Note that 
without loss of generality the divergent behavior was as- 
cribed to the first eigenvalue. However, it is clear that 
any particular enumeration of A M has no invariant mean- 
ing. What is actually meant in Eq. i|ll|) is that only one 
eigenvalue is quadratically divergent, but it is not possi- 
ble to assign a particular number to it before averaging. 
In order to convince the reader that (|llfl is compatible 
with both @ and (fTOf) we note that since Ao and A, de- 
pend upon completely different scales it is legitimate to 
write £ (A M A„) = (A ) £ (Ai) + £ (A»Aj). Then Eq. @ 
becomes 

a(3 -a^ 2 + ^2(X i X ] ) = A + B ■ a' 2 (12) 

and under quite natural (in view of (|llf> ) assumption 
£ (AjAj) ~ a qcd ^ leads to indeed stringent relation 
between various coefficients 

a-0 = B. (13) 



It is evident that this equality is highly non-trivial and is 
not guaranteed a priori. If confirmed by lattice measure- 
ments it would imply the validity of the ansatz (|11|) thus 
providing the mean to check its self-consistency. As is dis- 
cussed in section Ull AI Eq. (|13fl is strongly supported by 
the measurements and is fulfilled with rather amazing ac- 
curacy. Moreover, in that section we also consider along 
the same lines the triple correlator £ A1<t/<p (A M A^A P ), 
for which the analogous to relation holds and which 
provides the additional consistency check similar to i|13H . 
Note, however, that in the latter case the numerical un- 
certainties are larger and for this reason we do not con- 
sider the triple correlator here. 

To summarize, we found that the ansatz l|llfl repro- 
duces precisely the observed pattern of leading power 
divergences ©, JSJl, (|10fl and is fairly compatible with 
measured local characteristics of the topology defining 
map S 4 — > HP 1 . Since even the relation (|13fl is strongly 
supported by the data, we are confident that Eq. (|llf> re- 
flects correctly the leading UV dependence of the map 
S 4 — > HP 1 , which therefore turns out to be highly 
asymmetric on average. It is worth to note that the 
standard picture implies (A M ) cx Aq CD meaning that 
unit topological charge is gathered at large (of order 
A QCd) distances. The singular behavior of one eigen- 
value (Ao) cc 1/a 2 signifies immediately that the topolog- 
ical susceptibility \ = (Q 2 ) /V is saturated on submani- 
folds with characteristic four- volume of order a ■ ^qcd- 
In other words the topological density is to be concen- 
trated mostly in three-dimensional domains embedded 
into Euclidean four-dimensional space (we return to this 
problem in section III Bl) . 

The above results have rather dramatic consequences 
for the topological density q x , to illustrate which let us 
consider (q 2 ) defined as (q 2 ) = lim| 2; |_ > o(9o92;)- Note that 
within the usual approaches this definition is, in fact, 
ambiguous since the correlation function (qoq x ) is pertur- 
batively dominated at small distances (qoq x ) ~ — 1/ |a; | 8 . 
Hence the perturbative ambiguities make the value of 
(q 2 ) undefined and the standard lore is to fix it by the 
requirement \ = JiQoQx)- However, it is crucial that 
the HP 1 embedding approach to the gauge fields topol- 
ogy is factually exempt from perturbation theory. The 
best illustration is provided by Eq. J2J and in connection 
with topological density correlation function we discuss 
this in section Ull CI Therefore for HP 1 -based topolog- 
ical density, which in the continuum limit is given by 

q cx TrF HP F HP , the estimate (q q x ) l/|a;| 8 is not 

valid and UV behavior of (q 2 ) should be considered anew. 
Note also that the possible UV divergence of (q 2 ) could 
not be subtracted as usual and is not equivalent to con- 
ventional contact terms. From now on and throughout 
the paper the topological density q x is always understood 
via HP 1 embedding approach. 

In fact, the leading UV behavior of HP 1 -based topo- 
logical density could be established from Eq. J7J. Indeed, 
from the well known properties of instanton solution it 
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follows that 



q oc Tr F KP F HP K det J = ±det 1/2 



9- 



(14) 



Here it is convenient to introduce the notion of char- 
acteristic topological density q, which could be defined 
rigorously as q 2 = (q 2 )- However, below we'll sometimes 
use this term and the same symbol q to denote just the 
typical scale of the topological density fluctuations. The 
justification is that as far as the dependence on UV and 
IR scales is concerned these definitions essentially coin- 
cide. Therefore, from Eq. (|14fl we conclude that 



= (4 , )<x(dBtg) = ([[\ lJ 



A 6 

QCD 



(15) 



where Eq. Ijlljl had been used. The conclusion is that 
even the non-perturbatively defined topological density 
is divergent in the continuum limit reflecting directly the 
highly asymmetric local structure of the topology defin- 
ing map S 4 — > HP 1 . Note however that this divergence 
is incomparable with usual perturbative one (q 2 ) ~ 1/a 8 , 
which, in fact, is even non-integrable. Of course, it is un- 
derstood that the very definition of the topological den- 
sity is arbitrary to large extent (full derivative could al- 
ways be added). However, the HP 1 embedding method 
is specified completely with no free parameters involved. 
Moreover, the corresponding topological density is defi- 
nitely exempt from perturbative ambiguities so that the 
divergence l|15[) could not be equivalent to the contact 
term and should be dealt with accurately. In particular, 
we argue in the next section that Eq. I|15fl implies rather 
peculiar geometry of vacuum topological fluctuations. 



B. Dimensionality of Topological Fluctuations 

The problem to be addressed in this section is the ge- 
ometrical properties of vacuum topological fluctuations 
and in particular their dimensionality. First, let us note 
that any particular distribution of topological density in 
finite volume V could be divided unambiguously into the 
regions V+ (q x > 0) and V- (q x < 0) of sign-coherent 
topological charge so that the total charge is given by 



Q= q = 



- Q- 



(16) 



where the relative sign of Q + , Q_ was made explicit in 
the last equality (Q_ > 0) and no specific properties 
of the regions V+, V- were supposed (V± could be dis- 
connected, in particular). Note that in this section we 
do assume that the charges Q± could be estimated as 
Q± = qV±. For mean squared topological charge the 
decomposition H16|l implies 



(Q 2 ) « (Q\) 



(Q+Q-) 



(17) 



where we have generically assumed that {Q 2 _) — (Q+). 
Since the topological susceptibility x is « priori postu- 
lated to be finite in the continuum limit, we have 



(Q 2 + )-(Q + Q.) _ 2 (V 2 )-(V+V-) 



V 



V 



A 4 

1V QCD 



(18) 



Keeping in mind the ultraviolet divergence of q, Eq. 115fl . 
we conclude that the independent fluctuations of the vol- 
umes V± are strictly prohibited. Instead the magnitudes 
of V + and V_ are to be tuned up to the order 0(a 2 ) 



((V+-V-) 2 ) 



i 2 k- 2 



QCD 



■v, 



(19) 



since otherwise the topological susceptibility would di- 
verge in the continuum limit. It is important that, con- 
trary to the case of zero-point fluctuations, for non- 
perturbatively defined topological density there are no 
arguments which would guarantee the exact generic can- 
cellation of divergent terms in the integral Q = j q. In 
fact, the fine tuning assumption i|19[) is not new, simi- 
lar in spirit observations were made already in the recent 
past (see, e.g., Refs. [1, H, S E|)- However, we would 
like to reformulate the problem so that explicit powers of 
lattice spacing do not appear. Indeed, Eq. (|19|) implic- 
itly assumes that the volumes V± are four-dimensional 

and is satisfied identically if the topological density is 

( 

distributed in three-dimensional domains V± = V±/a, 
which are allowed to fluctuate on the scale of Aqc d 



((V 



(3) 



yW)2) „ v/A 



2 

QCD 



(20) 



Note that l|20|) is not the real solution, but rather the 
reformulation, of the fine tuning problem. Indeed, al- 
though the explicit spacing dependence is gone, the three- 
dimensional structure of topological fluctuations in D=4 
YM theory is equivalent, in fact, to a sort of fine tuning. 

In order to make the presentation more coherent, let 
us return to Eq. I|17fl . We could easily bypass Eq. i|19|) 
assuming that the charges Q± fluctuate independently 



(Q+Q-) = (Q + )(Q-) = (Q+y 

Then Eq. I|18|l translates into 



(Ql) - (Q^ 



V V " V 

where the validity of central limit estimate 



A 4 

1V QCD 



(Q' 2 



(Q+f cx (Q H 



(21) 

(22) 
(23) 



was supposed. Then the relation between ultraviolet be- 
havior of characteristic topological density and the di- 
mensionality of the corresponding fluctuations could be 
given as follows. Assume that q is of order a~ a and that 
the dimensionality of the topological fluctuations is D so 
that V + = a 4_£) ■ V+ , where vj ^ is spacing indepen- 
dent. Then 



,4-a-D 



{V. 



{D)\ 



V 



A 4 

1V QCD 



(24) 
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and the relation between UV behavior of the character- 
istic topological density and the dimensionality of the 
relevant topological fluctuations follows 

q ~ a- a , dim[V±] = 4- a. (25) 

However, it is clear that the argumentation relies heav- 
ily on the assumption that the fluctuations of the topo- 
logical charges Q± obey Eq. (|2*T|> (as well as Eq. 
which, however, seems to be less restrictive). A priori 
Eq. {2B is by no means evident and being confronted 
with experimental lattice data provides the most strin- 
gent test of the above scenario. Various experimental 
aspects of the problem are addressed in the next sec- 
tion. Here we note only that Eqs. 1)15(1. (|21J) imply the 
three-dimensional structure of vacuum topological fluc- 
tuations. Evidently, the reversed argumentation could 
also be given, namely, the dimensionality of sign-coherent 
topological charge fluctuations determines the leading 
ultraviolet behavior of characteristic topological density 
provided that Eq. I(21() is valid. We stress that the essence 
of the above presentation is the factual absence of lead- 
ing perturbative divergences in HP 1 projected fields and 
in the corresponding topological density. One could con- 
vince oneself that in the case of perturbatively dominated 
topological density, q ~ a~ 4 , the assumption 1(19(1 (with 
a 2 replaced by a 4 on the r.h.s.) holds true, while Eq. 1)21(1 
is violated and reads instead (Q+Q-) = (Q±)- The con- 
clusion is that Eqs. ((151) . ((211) . 1(231) are indeed crucial 
to validate the lower dimensionality of the topological 
charge fluctuations. 

III. CONFRONTING WITH LATTICE DATA 

In this section we describe in details the results of our 
numerical investigations of the scenario outlined above. 
In section IlII Al we study the lattice spacing dependence 
of the eigenvalues of the metric associated with the topol- 
ogy defining map S 4 — ► HP 1 and discuss in details the 
results announced in section III Al In section IIII Bl the 
topological density at various scales is considered; we 
show that even the simplest approach indeed qualita- 
tively confirms the divergence of characteristic topolog- 
ical density in the continuum limit. Section Till ("l is de- 
voted to the investigation of (qoq x ) correlation function 
from which we deduce the scaling law of the character- 
istic topological density. Then in section IIII Dl we pro- 
pose a method, which includes essentially no free parame- 
ters and allows to directly establish the dimensionality of 
topological fluctuations. Finally in section llll El the topo- 
logical charges associated with sign-coherent domains are 
shown to fluctuate indeed independently thus providing 
the self-consistency check of the above scenario. 

The numerical measurements were performed on 7 sets 
(TableQJ of statistically independent SU(2) gauge config- 
urations generated with standard Wilson action. The 
most of configurations listed in Table [I] are the same 
as were used in Refs. 0, |n| (except for the set at 
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a,fm 


L t 


L s 


V phys ,ha 4 


jqconf 


Ajconf 
iv <J 


2.4000 


0.1193(9) 


16 


16 


13.3(4) 


198 


70 


2.4273 


0.1083(15) 


16 


12 


3.8(2) 


250 


80 


2.4500 


0.0996(22) 


14 


14 


3.8(2) 


200 


80 


2.4750 


0.0913(6) 


16 


16 


4.6(1) 


380 


75 


2.5000 


0.0854(4) 


18 


16 


3.92(7) 


200 


75 


2.5550 


0.0704(9) 


20 


20 


3.9(2) 


80 


80 


2.6000 


0.0601(3) 


28 


28 


8.0(2) 


65 


60 



TABLE I: Simulation parameters. 



(3 = 2.555). The last column in Table |U represents the 
number of configurations on which we calculated the bulk 
topological charge density. Note that the number of an- 
alyzed configurations at each spacing is indeed rather 
large, which is due to the new algorithm used to eval- 
uate the topological density (the algorithm is described 
in Appendix). The lattice spacing values quoted in Ta- 
ble[I]were partially taken from Refs. [l|j and fixed by the 
physical value of SU(2) string tension ^/a = 440 MeV. 
Note that for (3 = 2.4273 and f3 = 2.555 the lattice 
spacings and corresponding rather conservative error es- 
timates were obtained via interpolation in between the 
points quoted in [Tfij ]. 

A. Local structure of the map S 4 — > HP 1 

The local structure of S 4 — > HP 1 mapping is charac- 
terized by the corresponding Jacobian Jjf (see ||5J, (0), 
1(14(0 . however in this section we concentrate on the in- 
duced metric g^ v = J-fJ^ and its strictly positive eigen- 
values X/j,, ji = 0, 3. As we noted already, any par- 
ticular enumeration of A M has no invariant meaning so 
that the meaningful observables associated with g^ v can 
not depend upon the ordering. Let us first describe the 
actual numerical procedure used to extract the spectrum 
{A^} at any particular lattice point x, at which we know 
the unit five-dimensional vector € HP 1 , A = 0, ...,4 
as well as the analogous quantities at the neighboring 
sites n£ +v S HP 1 . In accord with what had been said 
in section fll Al we introduce stereographically projected 
coordinates yg, y x+v such that y£ = and then consider 
the discretized approximation to the Jacobian 

Jf}{x) = o" 1 •(!£+„ - y£) = a- l -y» x+v , (26) 

from which the metric and its spectrum {A M } are 
obtained straightforwardly. It goes without saying 
that we confronted the quantity E M <„ X^\ v ](x) with 
[ \ Ti(Fj^') 2 ](x) on each configuration and found that 
Eq. JSJ is satisfied almost identically at every lattice site. 

The lattice spacing dependence was measured for three 
ordering insensitive observables associated with g^ v 

Mi(a) (27) 
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a 2 ,fm 2 
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0.012 



FIG. 1: Scaling of the correlators M n , Eqs. fM . 
with diminishing lattice spacing. Lines represent the best fits 
according to Eq. 13UII . 



M 2 (a) = ( x ^) 



Ms(a)= ^ X - X r- 



(28) 



(29) 



ii<v<p 



The results of our measurements are presented on Fig. ^ 
and indicate strongly that UV behavior of all these quan- 
tities is well described by 



M r 



A 2 ™ 

QCD 



Pn~ A 



1,2,3, 

2(n-l) 
QCD ■ 



(30) 



It is remarkable that the observed pattern of power cor- 
rections is universal and for each correlator (|27|l . (|28H . 
(12131) includes only terms of order O(a ) and 0{a~ 2 ). 

What is even more important here is that the numeri- 
cal data l|30|) is in perfect qualitative agreement with the 
assumption Moreover, the actual numerical val- 

ues of the coefficients a n , (3 n provide the self-consistency 
check of the ansatz 111 H . Indeed, considering Eq. (|30l) for 
n = 1 and n = 2 we obtain 



Oil ■ &1 = (h, 

2j(AiAy) = a 2 , 



(31) 



while the combination of n — 2 and n — 3 cases leads to 
& = Pi-J2(^ X j) = A- a 2. (32) 



Evidently, Eqs. (JSU, (E2J) are highly non-trivial and their 
numerical validity would definitely signify the correctness 
of the ansatz . As far as the numerical values of the 



relevant coefficients are concerned, 
best fits according to Eq. 130|) is 



the outcome of the 



a 2 ■ (X\ 
4 

a 1 ■ a 2 



11.8(3), 
87(3) , 



/3i = 0.018(2), 
i 2 -(3 2 = 0.24(2), (33) 



2.1(2). 



One can see that these numbers are in the perfect agree- 
ment with Eq. (|3*TJ> 



Qi • Pi 

P2 



0.9(1), 



and are compatible with Eq. (|32() 



0.8(2). 



(34) 



(35) 



although the deviation from unity and numerical uncer- 
tainty is larger in the latter case. 

The conclusion is that numerical data supports 
strongly the assumption l|ll(l so that indeed the topology 
defining map S 4 — > HP 1 is likely to be highly asymmet- 
ric. One of the relevant eigenvalues is divergent in the 
continuum limit and seems to depend only on the ultravi- 
olet cutoff. On the other hand, the remaining eigenvalues 
are sensitive solely to the infrared scale and show no di- 
vergences near the continuum limit. 

Finally, we note that it is tempting to consider along 
the same lines the fourth order correlator M4 — ([\ A M ) 
and then relate via Eq. (jl5[) its spacing dependence with 
leading UV behavior of the topological density. In turns 
out that numerically Ai^ia) is indeed compatible with 
(|30|l . however, we refrain to rely on (|15fl at finite lattice 
spacing The point is that Eq. i|15|) is certainly valid in the 
limit a — > 0, provided that Mi stays constant in physical 
units. However, due to the suspected power divergence 
the corrections to Eq. (|15fl are difficult to estimate. 



B. Topological Fluctuations at Various Scales 

As was repeatedly stressed in Refs. any discus- 

sion of the topological charge density within the lattice 
settings inevitably introduces a particular cutoff A q on 
the magnitude of the density so that q(x) is equated to 
zero if \q(x)\ < A q . Indeed, the most straightforward ar- 
gument here is that in the numerical simulations the den- 
sity is always known with finite accuracy. Thus the nu- 
merical precision provides the finest possible cutoff which 
in physical units evidently scales like A q oc a -4 . More- 
over, the introduction of (often implicit) finite A q is in- 
herent to all studies of the gauge fields topology. For 
instance, the chiral fcrmions based topological density, 
which is given by the sum of Dirac eigenmodes con- 
tributions, is usually either restricted to lowest modes, 
A < A oc Aqcd, or is considered for all modes available 
on the lattice, A < 1/a. Therefore the actual problem is 
not the presence of the cut A g , it is introduced always. 
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FIG. 2: Volume density of points with \q x \ > A 9 as a function 
of lattice spacing at different A 9 scaling laws 137H . 



The physically meaningful question is the spacing depen- 
dence A q = Aq(a) and the above examples illustrate two 
extreme cases A q cx a~ 4 and A g cx Aq CD . 

It is crucial that the scaling law A q (a) could be taken at 
will and we're going to exploit this freedom to study the 
spacing dependence of the topological density. Indeed, if 
the characteristic topological density q stays constant in 
physical units then the volume density of points at which 

\q(x)\ > Ag, 




\q(x)\ > A q 
otherwise 



(36) 



should also be lattice spacing independent for A g oc 
Aq Cd . Note that the lattice units had been used in 
(|36[) and that p{A q ) is dimcnsionless, positive, bounded 
p(A q ) < 1 quantity, unrelated to the volumes V± dis- 
cussed in section Hi Bl It is clear that the divergence q ~ 
a~ a would result in the divergent behavior p(A q ) ~ a~ a 
of the volume density |l9j for physical cut A q oc Aq CD , 
while p(Aq) is to be almost spacing independent for sim- 
ilarly divergent cut A g oc Aq CD ■ (a Aqcd) 01 ■ There- 
fore the most straightforward way to analyze the spacing 
dependence of the characteristic topological density is to 
tune the scaling law A q (a) until the volume density p{A q ) 
becomes constant at various lattice resolutions. 

Unfortunately, this approach does not allow to inves- 
tigate the dependence q(a) precisely. Indeed, on the lat- 
tice we could only probe a finite set of scaling laws A q (a) , 
moreover the corresponding estimates of p(A q ) are always 
biased. However, it is crucial that the correct qualitative 
picture could easily be obtained this way. We performed 
the measurements of the volume density p(A q ) at various 
lattice spacings using three different scaling laws 



A, n) cx A% CD -{aA QCD )- n . 



0,1,2, 



(37) 



where the numerical coefficients were chosen in such a 



(this particular choice is motivated below). The re- 
sults of our measurements are presented on Fig. [21 from 
which it is clear that the volume density of points with 
| q x | > Aq '' = [200 MeV] 4 is rapidly rising with dimin- 
ishing lattice spacing. Contrary to that the linearly di- 
vergent cut on the topological density, Aljp ~ l/a re- 
sults in the almost spacing independent volume density 
p(A q 1] ) w 0.35. On the other hand, once the quadrati- 
cally divergent cut A q 2 ^ ~ l/a 2 is imposed the quantity 
p(A q ) diminishes almost linearly with vanishing lattice 
spacing and in the limit a — > becomes compatible with 
zero. The conclusion is that the characteristic magnitude 
of the topological density is indeed singular in the con- 
tinuum limit, the leading divergence is compatible with 
linear one 



q oc 



Aqcd / a 



(38) 



way that A q n ^ 



[200 MeV] " 



0, 1,2 at a = 0.1193 fm 



and is in accord with theoretical expectations l|15|) . 

Let us discuss now the dimensionality of topological 
charge sign-coherent regions. Qualitatively the structure 
of topological fluctuations at various cuts A g is rather 
simple [ToL 111] ]. For utterly small values of A g there 
are typically only two large (percolating) regions of sign- 
coherent topological density, each of which occupies al- 
most half of the lattice volume and carries rather large 
topological charge Q±, Eq. (|16|l . With rising cutoff the 
volume density of percolating regions diminishes, while 
the number of small sign-coherent lumps rapidly grows. 
Finally a sort of percolation transition happens at which 
the percolating regions cease to exist and become indis- 
tinguishable from the small lumps. After that point the 
volume distribution of sign-coherent regions become uni- 
versal (Aq independent) and is described by rather re- 
markable power law. However, the physics changes dras- 
tically at the lumps percolation transition. Namely, the 
string tension, associated with HP 1 projected fields and 
which accounts for the full SU(2) string tension in the 
continuum limit, vanishes. Already from this observation 
we expect that the most physically important topological 
fluctuations are represented by the largest (at given cut- 
off A q ) lumps in topological density and it is natural to 
focus exclusively on their dimensionality. Note that the 
actual values of the cuts l|3*7|l were taken to be always be- 
low the lumps percolation transition in the whole range 
of lattice spacings considered. 

However, at any fixed cutoff A q the structure of the 
lumps is very complicated and their dimension is, in fact, 
not a well defined concept. Indeed, the notion of dimen- 
sionality makes sense only as a scaling relation since at 
any fixed lattice spacing the lumps occupy some finite 
fraction of the volume. Actually the situation is much 
worse since the very definition of the lumps require in- 
troduction of the cutoff A g , the spacing dependence of 
which is not fixed. Moreover, admitting the lower dimen- 
sionality of sign-coherent regions, their volume fraction 
is not obliged to be finite in the continuum limit, hence 
even the spacing dependence of lumps localization degree 
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(which might be expressed in terms of inverse participa- 
tion ratio or similar quantities) would not reveal their 
dimensionality. 

It is clear that the crucial obstacle in lumps dimension- 
ality definition is the necessity to impose the cutoff on the 
topological density. The concept of the dimensionality of 
topological fluctuations would become unambiguous pro- 
vided that we could get rid of explicit A g and hence reject 
the language of the lumps. This program is implemented 
in section UlI Dl However before going into details let us 
study the divergence 1)38(1 more quantitatively and con- 
sider the topological density correlation function. 



C. (qoq x ) Correlation Function 

It is was discussed in brief in section III Al that con- 
sidering the magnitude of the characteristic topological 
density defined by q 2 = lim^^o (qqIx) one has to prove 
that q indeed makes sense and is not equivalent to the 
contact term inherent to the perturbation theory. It was 
stated without proof that this is the case for HP^based 
topological density. In this section we present the cor- 
responding data and investigate the correlation function 
(loQx)- Note that the correlator (qoq x )js known to be 
negative at any non- vanishing distance [l4[ provided that 
the definition of the topological density is local (see also 
Refs. |f|,l3l|j| f° r discussions). However, the requirement 
of locality is a priori violated in HP 1 embedding approach 
so that (qaq x ), \x\ ^ is not obliged to be negative. We 
could only hope that the intrinsic non-locality is not so 
violent and extends up to some distance Rq fixed in phys- 
ical units. Note that this expectation is not completely 
groundless. Indeed, many non-perturbative observables, 
defined via HP 1 projection and studied in Refs. [Tol fTT| . 
do not reveal any pathology and reproduce, in fact, the 
corresponding results in the full theory. Actually the de- 
gree of non- locality of HP 1 method could be estimated 
from the behavior of heavy quark potential measured at 
a = 0.0601 fm in and it turns out to be of order 
Ro < 0.2 fm. On the other hand, the investigation of 
(lolx) correlation function allows to find Rq rather pre- 
cisely and check its scaling properties. 

Generically we expect that the correlator (qoq x ) is to 
be positive up to the distance Rq and then should be- 
come negative provided that Rq is finite. Remarkably 
enough these expectations are precisely confirmed by the 
measurements. The positive core of (qoq x ) correlation 
function at small is presented in physical units on 
Fig. [3J It is apparent that the points at various spacings 
are well described by the exponential |a;|-dependence 



-\x\/R qq 



<R a 



(39) 



Note that due to the logarithmic scale used on this plot 
the data sets are terminating at the same physical dis- 
tance, for larger \x\ the correlation function becomes neg- 
ative. Therefore the degree of non-locality inherent to 
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HP 1 embedding method is given roughly by 

Ra < 0.4 fm (40) 

and is indeed constant in physical units as is evident 
from Fig. 01 Moreover, the almost perfect scaling of var- 
ious data sets indicates once again that HP 1 projected 
fields do not contain any trace of the perturbation the- 
ory. Indeed, the mixture with perturbative contributions 
would lead to notable ~ — I/|a;| 8 terms and would result 
in rather abrupt deviation from the exponential behav- 
ior. On the other hand, at large distances, |x| > Rq, the 
correlation function (qoq x ) indeed becomes negative as is 
illustrated on Fig. 0] It is important that the negative 
part of (qoq x ) correlator does not show any singularity in 
the limit a — >■ 0, in particular, it has nothing to do with 
usual perturbative dependence. 

Let us consider the scaling properties of the correlation 
length R qq , Eq. (|39|l . As might be expected already from 
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0.12 



FIG. 5: Correlation length R gg , Eq. I39II . as a function of 
lattice spacing. Line represents the best linear fit. 



our measurements are presented on Fig. [B] from which 
it is clear that HP 1 -based characteristic topological den- 
sity still diverges in the continuum limit. However, this 
divergence has nothing in common with perturbatively 
expected 0(a~ 8 ) behavior, in fact, it is much weaker and 
is compatible only with linear or quadratic dependence 



(I 2 ) 



b„ 



1,2. 



(43) 



Unfortunately, our data points do not distinguish these 
power laws and are adequately described by either linear 
(Xn=i = 0-9) or quadratic {Xn=2 — 1-5) °ne. Neverthe- 
less, the ultraviolet divergence of characteristic topolog- 
ical density q ~ a~ a within the HP 1 embedding method 
could be considered as firmly established. Moreover, we 
are confident that the corresponding power exponent is 
close to unity, a < 1, and is in accord with theoretical 
expectations (tTsll and the estimate (|3*5|) obtained earlier. 
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FIG. 6: The mean squared topological density (q 2 ) as a func- 
tion of lattice spacing. Curves represent the best power law 
fits, Eq. |S3) . 

Fig. [3J R qq decreases with diminishing lattice spacing. 
In more details the dependence R qq (a) is presented on 
Fig. [5J from which it is apparent that the correlation 
length is likely to be linear function of a with rather small 
continuum value 

R qq = 0.037(1) fm. (41) 

As is apparent from the above presentation, the short 
distance behavior of (qoq x ) correlation function is defi- 
nitely exempt from the perturbative uncertainties and in 
the limit |x| — > the characteristic topological density 

f = (q 2 ) = lim (q q x ) (42) 

has nothing to do with contact terms inherent to the 
usual approaches. Therefore let us consider the depen- 
dence of q upon the ultraviolet cutoff. The results of 



D. Dimensionality of Topological Charge 
Fluctuations 

1. Biased Random Walk Model and the Choice of 
Parameters 

In this section we introduce the model, which allows 
to investigate directly the dimensionality of the relevant 
topological fluctuations. Although the resulting method 
is not entirely rigorous, its ambiguity reduces to only one 
free parameter, the choice of which we thoroughly dis- 
cuss. Generically the idea is to consider some dynamical 
system the evolution of which is sensitive to the dimen- 
sionality of the ambient space. Then if we embed some- 
how this system into the topological density background, 
its evolution will reveal the effective number of available 
dimensions, which is to be naturally associated with the 
dimensionality of the relevant topological fluctuations. 
The simplest dynamical system of this sort could be con- 
structed on the top of usual diffusion equation, which in 
turn is equivalent to the model of random walks. Then 
the dependence upon the external environment could be 
introduced by making the hopping probabilities to be the 
local functions of the ambient space characteristics; the 
models of this sort are known as biased random walks. 
Correspondingly, the effective dimensionality D as is seen 
by biased random walkers is referred to as diffusion di- 
mension. The purpose of this section is to precisely for- 
mulate and investigate this approach. 

We start from the observation that the introduction of 
the cutoff A q on the topological density (see section llll Bp 
was aimed solely to separate the inevitably present noise 
(utterly small values of q x ) from the relevant fluctua- 
tions, which are associated with relatively large values of 
q x . Although the choice of particular A q indeed makes 
the notion of 'small' and 'large' well defined, the geome- 
try of the resulting lumps in topological density strongly 
depends upon the cut. It is apparent that the weak point 
here is sharpness of the dimension four cut, it would 
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be much advantageous to remove the small topological 
density regions softly, making the lumps geometry much 
more robust with respect to the parameters involved. It 
turns out that slightly modified diffusion model is indeed 
suitable to achieve this. Namely, we propose to mod- 
ify the diffusion equation by allowing the random walk 
to hop towards the regions of higher topological density 
with larger probability. In the language of the diffusion 
equation, which describes the propagation of heat, this 
amounts to the introduction of space-time dependent dif- 
fusion coefficient (thermal conductivity), which vanishes 
in the regions of small topological density, and hence heat 
is allowed to spread only within the domains of large \q x \. 
Then the decay rate of the initial heat pulse, which is the 
same as the return probability for corresponding biased 
random walk, essentially reflects the number of available 
dimensions within the topological fluctuations and hence 
is to be identified with their dimensionality. 

In fact, this general idea fixes almost uniquely the ran- 
dom walk model which we would like to investigate. It is 
convenient to start directly from the microscopic rules of 
the biased random walk, which require that the proba- 
bility p Xy x+fi to hop from point x to the neighboring site 
x + fi is monotonically rising scale free function of the 
topological density magnitude at x + fi 

V , = \^+^ (AA) 

where the power exponent 7 > could not be fixed a pri- 
ori and remains free parameter of the model. Note that in 
this section we exclusively consider the local magnitude 
\q x \ of the topological density; to lighten the notations 
the corresponding modulus sign will be omitted. From 
Eq. Ij44() it is straightforward to obtain the continuum 
diffusion-like equation, which determines the probability 
P(t, x) to reach the point x during the proper time inter- 
val t provided that at t = walker starts at x — (see, 
e.g., Ref. 0) 

d t $ = \q- 2 ~'d x [q 2 "<d x <S>] , (45) 
^(t,x) = P(t,x)/q^(x), 

where the initial condition is P(0,x) = 8(x). However, 
Eq. H45f) is not yet the usual diffusion equation, in par- 
ticular, the decay rate of the initial perturbation (return 
probability in the random walk language) is not guaran- 
teed to be P(t, 0) oc t~ D / 2 . To get the correct interpreta- 
tion of H45f) we introduce new coordinates Q{x) according 
to 

d^jdx v = q^-S^. (46) 

In the particular case of one dimensional problem l|4t)|) 
allows explicit solution C,(x) = f x q 2j , where the lower 
integration limit is taken at arbitrary fixed point. It 
is important that £ is a single valued function of x al- 
most everywhere, moreover, its range is determined by 
the magnitude of the topological density. Indeed, the 



regions of utterly small \q x \ are squeezed to almost one 
point by the map H46|) regardless how large these regions 
were in x space. It is crucial that the term 'small' above 
obtains unambiguous and physically correct meaning of 
relative smallness since only the relative variation of the 
topological density does matter. Indeed, Eq. I|45(l is ev- 
idently scale invariant under q —* Xq and is equivalent 
to usual diffusion equation for everywhere constant q x . 
In terms of new coordinates Eq. (|45(l takes the standard 
form 

dt® = ld c [g 47 9 c *] ■ (47) 

o 

We conclude therefore that the diffusion process i|44|l . 
Q45J1 . (|47|l takes place in the regions of relatively large 
topological density and hence should reflect properly 
the dimensionality of underlying topological background. 
Moreover, the decay rate of the initial perturbation is 
given by 

$(t,Q) oc r D ^' 2 , (48) 

where D(^) is the diffusion dimension of the topological 
fluctuations and we have explicitly indicated that the di- 
mension defined this way depends non-trivially upon yet 
not fixed parameter 7 to be discussed next. 

The non-triviality of the dependence D(-f) is evident 
since in the limit 7 — > the model (|44|l . (|45|l reduces to 
standard unbiased random walk with 

£( 7 ^0) = 4, (49) 

while at 7 — > 00 the microscopic probability l|44|) allows 
the hopping only towards largest neighboring \q x \. Hence 
the walker is trapped eventually at the local maxima of 
topological density distribution and 

D( 7 ^oo) = 0, (50) 

regardless of the background. Already from this obser- 
vation it is apparent that without an additional physical 
input the model (|44|l . (|45|l would be essentially useless 
since at various 7 it reflects, in fact, different properties 
of the underlying background. For instance, the limiting 
behavior l|49|) implies that the topological density distri- 
bution is such that q x = only on measure zero set, 
while l|50f) is valid generically, provided q x is not identi- 
cally constant. Therefore the actual problem is not the 
arbitrariness of 7 parameter, but rather yet not posed 
physical question we're trying to investigate. 

In order to gain a physical insight it is crucial to re- 
tain the qualitative picture of vacuum topological excita- 
tions outlined in section ITlI Bl It was argued that phys- 
ically most important fluctuations are associated with 
percolating sign-coherent regions, the ultimate qualita- 
tive properties of which are the significant internal topo- 
logical density and extremely large linear extent. It goes 
without saying that we are interested precisely in the 
geometry of these sign-coherent domains. However, the 



11 



percolating lumps are not equivalent to just the regions 
of largest topological density as is revealed by the per- 
colation transition eventually appearing with rising A q 
cutoff. Even at large A q one finds individual 'hot spots' 
of small but non-vanishing volume which indeed possess 
the largest topological density. Consider now the decay 
rate of the initial heat pulse in the model (|45|l at large 
but finite 7, so that the effective thermal conductivity is 
non-zero only within the topological density hot spots. 
Evidently, in this regime the equilibration time t eq is fi- 
nite and is dictated by the typical size of the hot spots 
being much smaller than the squared lattice size L 



teq < £ 



for 7 > 7* 



(51) 



(the definition of 7* will become clear in the moment). 
Note that the distinct feature of this regime is that the 
double logarithmic plot of <I>(i, 0) significantly bends up- 
wards at t eq and therefore we expect generically that 



7 > 7 



d 2 ln$(i,0) 
d(\nt) 2 



>0, 



(52) 



reflecting, in particular, the drop in the effective thermal 
conductivity outside the hot spots. Evidently, at these 7 
values we are not probing the percolating lumps geome- 
try, the random walks are confined within the regions of 
largest topological density. 

Let us now gradually diminish the 7 parameter. The 
positive jump in the logarithmic derivative <91n<I>/<91nf, 
the location of which we still denote by t eq , would become 
smaller respecting the diminishing difference of thermal 
conductivities inside and outside the hot spots. Note 
that neither of the corresponding diffusion dimensions 
at t ^ t eq could be identified with the dimensionality of 
sign-coherent regions since the random walks are still too 
sensitive to the local irregularities of topological density 
even within the coherent domains. It might happen that 
at the particular value 7 = 7* the double logarithmic 
plot of Q(t, 0) degenerates into the straight line so that 



t eq w L 2 and 



7 = 7 



d 2 ln$(i,0) 
d(\nt) 2 



, t < L 2 



(53) 



Note that at this point the only relevant dimensional pa- 
rameter is the lattice size. Therefore the model becomes 
essentially scale free and its dynamics in the vicinity of 7* 
is reminiscent to the disorder driven conductor-insulator 
transition of condensed matter systems. It is crucial that 
at the critical point 7 = 7* the heat starts to propagate 
with constant rate over largest available distances, but 
since the thermal conductivity is significant only within 
the lumps, the heat transfer goes through the percolat- 
ing sing-coherent regions. In turn the condition 1(53(1 im- 
plies that local inhomogeneity of the topological density 
within the percolating lumps is inessential. Therefore it 
seems that only at the critical point 1(53() we indeed could 
obtain a consistent reflection of the relevant topological 
fluctuations in the biased random walks model. 



The existence of at least one critical value, 7* = 0, 
follows from 1)49(1 . but it is trivial one and surely ex- 
ists for arbitrary background q x . However, it is cru- 
cial that the above qualitative picture of vacuum topo- 
logical fluctuations hints on the existence of non-trivial 
critical coupling 7* > 0, which could arise entirely due 
to the low-dimensional long-range order present in the 
topological density distribution. For suppose that the 
topological charge sign-coherent regions are indeed lower 
dimensional objects possessing relatively large uniform 
topological density and extending through all the vol- 
ume. Then at the particular 7 = 7* the effective thermal 
conductivity would be non-zero only within the perco- 
lating regions, the lower dimensionality of which forbids 
the appearance of an additional scale in the heat prop- 
agation problem. Thus at this point the random walks 
model indeed contains no dimensional parameter apart 
from the lattice size and Eq. ((48(1 is fulfilled, while the 
corresponding diffusion dimension is to be identified nat- 
urally with the dimension of sign-coherent regions. Note 
that the lower dimensionality is crucial here, for regions 
of finite thickness the critical point 1(53(1 docs not exists. 
In fact, the analogous but reversed argumentation could 
also be given, namely, the appearance of non-trivial crit- 
ical point ((53(1 signifies the presence of lower dimensional 
long-range order in the topological density background. 

Note that the restriction to sign-coherent regions is au- 
tomatic in our approach. Indeed, as far as the model l|43|) 
is concerned, it was considered in the continuum limit as- 
suming differentiability of q x . Hence the domains q x ^ 
are separated by regions with vanishing thermal conduc- 
tivity. On the lattice the situation is more involved, but 
since the HP 1 -based topological density is definitely ex- 
empt from perturbative contributions we could be almost 
confident that the same argumentation applies. 

To summarize, the proposed biased random walk 
model seems to be efficient to investigate the structure of 
percolating topological density regions only at the criti- 
cal points, at which its dynamics becomes effectively scale 
free. We argued that the qualitative picture of the rel- 
evant topological fluctuations obtained earlier suggests 
the existence of non-trivial critical point, at which the 
diffusion dimension is to be identified naturally with the 
percolating regions dimensionality. 



2. Diffusion Dimension of Percolating Lumps 

Prior to presenting the results of our measurements 
let us discuss the lattice specific features of the above bi- 
ased random walk model, which complicate the numerical 
evaluation of the percolating lumps dimensionality. For 
any particular 7 value the solution of Eq. IpEST) could be 
constructed straightforwardly by implementing the ran- 
dom walk process, the rules of which are completely spec- 
ified by ((44(1 . The only subtlety here is the choice of 
the random walk starting point. Indeed, to improve the 



12 



0.01 



0.001 r 



© 0.0001 r 



1e-05 



1e-06 



Y= 0.10 i — h 
Y=0.25 
y 0.27 >■-* 
Y= 0.30 
Y= 0.80 





1e-05 r 
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Measurements were performed on j3 = 2.500 data set. 



statistics it is desirable to consider 

Mo $(o, Co; Co) 



*(t,0) 



Mo 



(54) 



= / dx Q $(0,x o ;t,xo) 



J I' 
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where $(0, xq', t, x) is the probability to reach the point x 
during the proper time t starting at xq. Eq. I|54[l means 
that the random walk starting point in x space is to be 
taken with probability oc q 8j . 

Therefore the crucial problem is how to find the rele- 
vant 7* value numerically and estimate the accuracy of 
its determination. For any fixed 7 we measured $(t, 0) 
by first taking some random starting point, chosen with 
oc q Sj probability, and considering in accord with 1441) 



FIG. 9: Return probability <&(i,0), Eqs. calculated 
on p = 2.500 data set with random permutation of topological 
density records on each configuration. 





2.4000 


2.4273 


2.4500 


2.4750 


2.5000 


2.5550 


2.6000 


7* 


0.27(3) 


0.27(3) 


0.27(4) 


0.27(4) 


0.26(2) 


0.27(3) 


0.28(3) 



TABLE II: Estimated 7* values for various data sets from 
Table □ 



the walk of total length 8 • 10 3 . Then the quantity 
$(i, x) = P(t, x)/q 2 (x) was constructed for each random 
walk and averaged with respect to ~ V/2 different start- 
ing points per configuration (V is the lattice volume). We 
checked that the statistics is large enough so that our re- 
sults do not change if more random walks are considered. 
Fig.|7|represents the double logarithmic plot of the quan- 
tity Q(t, 0) obtained on our /3 = 2.500 set at three close 7 
values, 7 = {0.25,0.27,0.30}, as well as the same quan- 
tity measured at 7 = 0.1 and 7 = 0.8. Note that for read- 
ability reasons different graphs are slightly shifted with 
respect to each other. As far as the data points at 7 = 0.1 
and 7 = 0.8 are concerned they definitely correspond to 
the regimes 7^7* and 7^7* respectively. Indeed, the 
first graph significantly bends downwards at t ~ 10-^20 in 
apparent disagreement with both Eq. (|52|) and Eq. H53|) . 
However, it is demonstrated below that negative second 
logarithmic derivative of $(i, 0) is a generic feature of the 
diffusion model in random environment. Note that at 
7 = 0.1 the system indeed equilibrates eventually at t eq 
of order few hundred, which is similar to squared size of 
the lattice used in this calculation. On the other hand, at 
7 = 0.8 the system starts to equilibrate at t eq « 20 <C L 2 
and the corresponding graph has positive second deriva- 
tive in accord with Eq. I|52|l . 

The inspection of the intermediate 7 values reveals that 
7 = 0.25 data points essentially lie on one single line, 
while the data set at 7 = 0.30 seems to deviate upwards 
from linear dependence. Note that the consideration of 
7 = 0.27 points are rather inconclusive, apparently the 
data just begins to bend aside the pure power law. There- 
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fore the relevant 7* value seems to be located around 
7* ps 0.25. We could estimate 7* more rigorously by con- 
sidering the quality of power law (|48|) fits at different 7. 
We fitted our data to Eq. gHl m tne range t < 100, the 
resulting x 2 /n.d.f. values are presented on Fig. [HI It is 
apparent from this figure that fits favor the critical value 
7* = 0.26(2), where rather conservative errors coming 
from Fig. [S] are quoted. 

In order to convince the reader that graphs on Fig. 
are indeed directly related to the underlying geometry of 
topological fluctuations, let us consider the same model 
(|44|l , l|45|) with the same parameters, but in the genuinely 
random external environment. Such a background, which 
is the 'nearest' to one on Fig.[7| could be obtained by ran- 
dom permutation of the topological density records on 
each configuration. We did this with our = 2.500 data 
set and then performed the identical measurements of 
$(t, 0) at 7 = {0.10, 0.27, 0.80}. The corresponding dou- 
ble logarithmic plots are presented on Fig. El and differ 
violently from that on Fig. [7] The most crucial obser- 
vation is that the second derivative for all graphs stays 
negative implying that the reasoning of section ITlI D H is 
inapplicable for randomly permuted topological density. 
Of course, this is not surprising since all the coherent 
structures are gone in the present case. The only limit 
in which the second logarithmic derivative of 0) van- 
ishes is given either by Eq. (|4*9|) . D = 4, or Eq. 1(50) 1. 
D = 0. We conclude therefore that without random per- 
mutations the non-trivial critical point 7* arises most 
likely due to the presence of long-range order in the topo- 
logical density. 

We performed the measurements of the relevant 7* val- 
ues for all our data sets listed in Tabled It goes without 
saying that indeed in each case the biased random walk 
return probability at 7 = 7* strictly obeys the power law 
(|48|l . The resulting estimations of corresponding 7* val- 
ues are summarized in Table [H] It is remarkable that 
the estimates of 7* appear to be spacing independent 
well within the quoted numerical uncertainties. It is true 
that the presented errors might be overestimated, how- 
ever, we are confident that 7* (a) dependence could not 
be revealed by the above method. Apparently this is 
due to the fact that 7 parameter is dimensionless so that 
the violent power dependence on the lattice spacing is 
unlikely to appear. 

Once the relevant 7* values are determined for every 
our data set it is straightforward to estimate the dimen- 
sionality of sign-coherent percolating regions of the topo- 
logical density. By construction the corresponding fits to 
Eq. (|48|l are practically perfect for every j3 so that we 
do not present full details of the fitting procedure. The 
point which should be discussed, however, is D(p/*) error 
estimation. In fact, the uncertainty coming from the fits 
to Eq. ((4*S)l is negligible compared to the ambiguity in 7* 
values. Therefore the errors in D(^*) were obtained by 
repeating the fits to Eq. I|48|) at minimal and maximal 7* 
within the corresponding error bands. Finally, we note 
that all fits were performed in the range t < 100. 
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a, fm 

FIG. 10: Diffusion dimensions D(j*), Eqs. 1051, at 

various lattice spacings. Line represents the fit to the constant 
behavior. 

The dimensionality of percolating sign-coherent re- 
gions of topological density measured with the above 
procedure at various lattice spacings is summarized on 
Fig. 1101 Remarkably enough the diffusion dimension 
D("f*) seems to be almost independent upon the lattice 
resolution well within the uncertainties. Its continuum 
value could be estimated from fit to constant behavior, 
which indicates in turn that D(j*) is definitely smaller 
than 4. Instead the numerical data suggests strongly 
that the diffusion dimension of the relevant topological 
fluctuations is 

D = £>(7*) = 3.07(3), (55) 

which seems to be in full agreement with their proposed 
three dimensional structure. Eq. (|55|l is in agreement 
with both theoretical expectations and experimental data 
on the topological density. What is still to be consid- 
ered is the internal consistency of the scenario outlined 
in our paper, which is expressed by Eqs. J2J, 1(2^)) ; this 
is the subject of the next section. Note, however, that 
the diffusion dimension is not the only definition of the 
dimensionality. It is still worth to confirm Eq. I|55|l with 
other methods. 



E. Consistency Check 

As was discussed in section III Bl the crucial equations 
which relate the divergence of characteristic topological 
density q and the dimensionality D of the relevant topo- 
logical fluctuations are 

(Q+Q-) = (Q+)(Q-) , (56) 



(Ql) - (Q±f cx (Q±) , (57) 
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FIG. 11: Fluctuations of the topological charges Q± at var- 
ious lattice spacings. Circles: spacing (in)dependence of the 
ratio A(Q±), Eq. 1581 : squares: the magnitude of Q± relative 
fluctuations characterized by B(Q±), Eq. J&BJ. 



where the charges Q± are to be calculated without any 
cutoff A g imposed. Moreover, the spacing independence 
of 156(1 . H57fl (if confirmed by the data) allows to put 
rather stringent restrictions on both q and D. As far as 
the lattice measurements are concerned, Fig. 1111 repre- 
sents the ratios 



A{Q±) 



(Q+KQ-) 



B(Q±) 



(QD - (Q±r 

(Q±) 



(circles) , (58) 
(squares) , 



as a functions of lattice spacing, where to improve the 
statistics the generic equalities {Q\) = (Q 2 -), (Q+) = 
(Q-) were used in calculation of B(Q±). As is evident 
from that figure Eq. 1(56(1 is satisfied identically in the 
whole range of considered spacings while the proportion- 
ality coefficient entering Eq. lj57|l does not depend upon 
the lattice resolution well within numerical errors. There- 
fore the validity of Eq. I(25|) is firmly established. Let us 
summarize the emerging qualitative picture of vacuum 
topological fluctuations which arises from the numerical 
data restricted by (|25*|) . 

It is apparent that the most confidential data is avail- 
able for the characteristic topological density. The theo- 
retical arguments based on the existence of the quadratic 
correction to the gluon condensate as well as the obtained 
numerical data suggest strongly that the topological den- 
sity is divergent at most linearly in the continuum limit 



< 1. 



(59) 



On the other hand, the dimensionality of topological fluc- 
tuations is still not firmly established and is plagued by 
theoretical uncertainties. Various estimations made both 



in this paper and in the literature 0, 0> El IS El suggest 
that it is of order three 



dim[V± 



3. 



(60) 



It is remarkable that Eqs. (|59|l . 160(1 overlap only at one 
point consistent with Eq. I|25() 



a 



1 . 



dim[V± 



3. 



(61) 



Given that Eqs. 1(56(1 . 1(57(1 are fulfilled with amazingly 
high accuracy we are forced to the conclusion that 1(61(1 
is the only values consistent with both theoretical con- 
siderations and numerical data. 



IV. CONCLUSIONS 

In this paper we further developed the SU(2) gauge 
fields topology investigation method, based on the em- 
bedding of HP 1 a- model into the given gauge back- 
ground [lOl jTlj . Our prime purpose was to exploit 
the remarkable properties of HP 1 projected fields found 
in , namely, the factual absence of leading perturba- 
tive divergences and simultaneous existence of non-trivial 
quadratic power correction to the gluon condensate. It 
is clear that these striking features of the projected fields 
are to be encoded into the local structure of the topol- 
ogy defining map S 4 — > HP 1 and hence should be re- 
flected in the corresponding topological density. The ex- 
tended analysis of leading power corrections performed 
in this paper allows to conclude that the topological 
density is likely to be linearly divergent in the contin- 
uum limit. Note that this divergence has nothing to 
do with short distance perturbative singularities and is 
much weaker. The divergence of the topological density 
by itself is almost academical problem since it is not di- 
rectly observable. However, combined with the require- 
ment of ultraviolet finiteness of topological susceptibility 
it leads to rather dramatic consequences for the geome- 
try of relevant topological fluctuations. Namely, we ar- 
gued that the topological charge is to be concentrated in 
three-dimensional submanifolds of four-dimensional Eu- 
clidean space. This is the only conclusion compatible 
with physical topological susceptibility, numerically es- 
tablished pattern of power corrections and which does 
not necessitates explicit fine tuning of the Yang-Mills the- 
ory at UV scale. Moreover, the fine tuning assumption is 
testable and if it does not happen, then one could derive 
rather stringent relation between the divergence of the 
topological density and the dimensionality of submani- 
folds, which support the most of the topological charge. 
Note, however, that the lower dimensionality of topo- 
logical fluctuations could also be considered as a sort 
of fine tuning, in which the explicit powers of UV cut- 
off are traded for unusual geometric properties. Qual- 
itatively our results are in accord with modern trends 
in the literature, which discuss the lower dimensional- 
ity of physically relevant vacuum fluctuations 0, lJ and, 
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in particular, of the topological density sign-coherent re- 

gions nasi sim. 

The actual experimental verification of this scenario 
turned out to be rather intricate both conceptually and 
technically; we believe that all these problems were ad- 
equately addressed in our paper. The technical achieve- 
ment is the development of fast and rather precise numer- 
ical algorithm of topological density evaluation, which 
allowed us to investigate the problem on the convincing 
statistical level. While the UV behavior of the topologi- 
cal density could be studied directly, the dimensionality 
of relevant topological fluctuations is much more involved 
problem, which consists essentially in physical interpre- 
tation of the term "relevant" above. We argued that 
the natural approach is to embed a dynamical system 
into the topological density background, the evolution of 
which is sensitive to the dimensionality of ambient space. 
The simplest system of this sort could be constructed on 
the top of usual random walk model and depends upon 
one dimensionless parameter. Then the phase structure 
of the system and location of its critical-like points is 
ultimately related to the long-range properties of the un- 
derlying background and, in particular, to the dimension- 
ality of sign-coherent topological regions. 

As far as the results of numerical experiments are 
concerned, our data show unambiguously that non- 
perturbatively defined characteristic topological density 
is divergent in the continuum limit. Moreover, we were 
able to obtain the upper bound on its leading spacing de- 
pendence. At the same time, the dimensionality D of the 
relevant topological fluctuations was shown to be decid- 
edly less than four and is compatible with D = 3. Here 
the assumed absence of the fine tuning becomes crucial 
and we showed that it indeed does not happen. Instead 
the topological charges associated with sign-coherent re- 
gions fluctuate independently. We conclude therefore 
that the only possibility to satisfy all the restrictions is 
to have linearly divergent topological density distributed 
in three-dimensional domains. 
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Appendix 

Here we describe in details the numerical algorithm 
used in this paper to calculate the topological charge 
density. The general formulation of the problem could 
be found in Refs. [ToL Hi). Essentially it reduces to the 
evaluation of the volume V(T) of 4-dimensional spher- 
ical tetrahedron T embedded into S 4 with vertices 



o 
E 
> 




FIG. 12: Cumulative distribution of V M c(T) and V det (T) (see 
text). 



given in terms of five i ■ 
A = 0, ...,4 vectors, (r?i) 2 



= 0, 
= 1. 



, 4 unit five-dimensional 



It is clear that for near vanishing V(T) the volume is 
given by V(T) = det Ai[nf]. Thus the volume of finite 
tetrahedron could be found by triangulating it into the 
set of small tetrahedra and summing up the infinitesi- 
mal volumes. In fact, our method works recursively until 
the determinant estimation of the volume of input tetra- 
hedron is larger than 10 -6 ; this way we indeed obtain 
the optimal performance. However, the determinant- 
based volume estimation is reliable and uniform only if 
the angles between all input vertices are small enough, 
e.g. (rii,rij) > 0. In this case it suffices to place new 
triangulation vertex at m oc fti and complete the re- 
cursion cycle. However, this procedure must be modi- 
fied if there are at least two input vertices i and j for 
which (rii,rij) < 0. Indeed, in this case the above re- 
cursions converge non-uniformly and eventually lead to 
almost degenerate tetrahedra with large volume but still 
almost vanishing determinant. The needed modification 
is to take the new vertex at m oc fti+ftj which guarantees 
that eventually we will get (rii,fij) > Vi,j. 

Let us note that for given accuracy of the determinant 
volume estimation for small tetrahedra (which is 1CP 6 
in our case) the volume of the original tetrahedron is 
evaluated, in fact, with finite bias which is due to the 
sphericity of every small tetrahedron. However for small 
enough volumes the sphericity could be accounted for by 
simple rescaling of the determinant estimation. To cali- 
brate the present algorithm we compared it with our pre- 
vious Monte-Carlo based method. To this end we gener- 
ated 5- 10 4 random spherical tetrahedra and applied both 
algorithms to each of them thus obtaining Monte-Carlo 
Vmc{T) and determinant-based Vdet{T) volume estima- 
tions. Fig. E| represents the cumulative distribution of 
Vmc(T), Vdet(T), which turns out to be astonishingly 
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narrow and is fairly compatible with linear dependence 

V MC {T) = J-V det (T), 7 = 1.1285(5), (A.l) 

where the optimal value of 7 coefficient results from the 
best linear fit. Note that the plot on Fig. ^] is restricted 
to V(T) < 5 • 10~ 3 , which is far beyond the maximal 
value of topological density even for our largest spac- 
ing; however, the linear dependence (|A.1J| remains valid 
even at larger V(T). Thus we are confident that the 
new triangulation method is definitely compatible with 
old Monte-Carlo approach in the relevant range of lattice 
spacings. 



Finally, we performed the same check as one described 
in Namely, we confronted the global topological 

charge, which could be found unambiguously for each our 
configuration with one calculated with present algorithm. 
It turns out that they agree in all cases with no excep- 
tions. Moreover, we found that the determinant-based al- 
gorithm gives even narrower distribution of Q fi oat around 
integer numbers compared to that of old Monte-Carlo 
based approach. We conclude therefore that the new 
method of the topological density calculation is superior 
to the old one both in accuracy and performance. 
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